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ABSTRACT 



The factors affecting the stability and the maneuvera- 
bility of a dynamically unstable ship were studied using 
the linear and non-linear equations of the motion of the 
ship. The DSL/360 Language was used to simulate the dynam- 
ics and to study both unstable motions and standard 
maneuvers . 
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TABLE OF SYMBOLS 



&eft= effective aspect ratio 
A f = area of movable fin 



A 



R 



area of rudder 



Cp = drag coefficient 
= lift coefficient 

0 = origin of reference axes fixed in the ship 

]£). = midship section 

x,y,z = system of reference axes fixed in the ship 

x ,y ,z = system of reference axes fixed in the earth 
o o o 

x G’ y G ,Z G = com P° nents distance of C.G. of the ship 



x^ = distance from 0 to center of pressure of fin along 
x-axis . 

ki.kgjk' = Lamb's coefficients (see ref. [l]) 
k^kg = coefficients of rudder deflection rate. 

I ,1 ,1^ = mass moments of inertia of a ship about x,y and 



m = actual plus added mass in the longitudinal direction 

X 

non-dimensional ized . 
m = mass of ship 



along x, y and z-axis, respectively. 



x ’ y ’ z 



z-axis, respectively. 



m^ = lateral added mass, non-dimensionalized . 



m. 



x 



rotational added mass, non-dimensionalized. 



<f> = roll angle 



9 = pitch angle 




8 



tes t<^ 



/3 = angle of attack in yaw on the hull (or drift angle). 

S - angular deflection of rudder 

U = resultant velocity of the fluid relative to a control 
surface or a fixed fin. 

p,q,r = components of resultant angular velocity of the 
ship about the x, y and z-axis, respectively 

( = 0, e,y>). 



u,v,w, = components of U along the x, y and z-axis, 
respectively . 

D = drag force 
F = total force 

L = lift force or length of ship (LBP) 

X,Y,Z = components of resultant total force acting at 
the origin 0 directed along x,y and z-axis, 
respectively . 

X , Y , Z = values of X,Y and Z at initial equilibrium 
o o o 

condition . 

H = angular momentum 

H = draft of ship 

= angular velocity 

= resultant total moments 

K,M,N = components of a resultant total moment acting on 

a ship about the x,y and z-axis, respectively. 

Y^ = partial derivative of Y with respect to v 

y. = " " v 

v 

Y = ..... r 

r 



Y- = 
r 

Yr = 



r 

s 
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N v = partial 


derivative of 


N with respect 


to V 


N ■ = 


1 1 


1 1 


V 


V 


N = 


!1 


1 1 


r 


r 


N- = 


II 


1 1 


r 


r 


Nj = 


1 1 


1 1 
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Note : 

1) Symbol of prime ( ' ) , except where otherwise stated 
indicates the values of an item non-dimensionalized . 

2) Symbol of dot ( * ) indicates differential value of 
an item with respect to time. 

3) Subscript f indicates the value of an item of fin. 

4) Subscript h indicates the value of an item of ship's 
hull. 

5) Subscript f^ indicates the value of an item of the 
large deadwood included rudder. 

6) Subscript f^ indicates the value of an item of the 
small deadwood or rudder only. 
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INTRODUCTION 



I . 

A ship is said to be dynamically stable on a straight 
course or in a turn of constant curvature if, when slightly 
disturbed from it's steady motion, it will soon resume that 
same motion along a slightly shifted path without any cor- 
recting control being applied. A ship which is dynamically 
unstable in straight line motion cannot maintain straight 
line motion when there is no control. However, with regard 
to maneuverability, a ship which is too stable will not 
turn as tightly as a somewhat less stable ship so that a 
highly stable ship may compromise maneuverability. There- 
fore, a ship should be designed for a moderate amount of 
stability so as to be able to go straight without compro- 
mising maneuverability too much. 

With these factors in mind then, the stability of five 
different ships was studied in Part III and maneuverability 
in Part IV. The first objective of this thesis was to 
investigate the stability of a ship with movable fin. The 
second objective was to investigate the maneuverability of 
an unstable ship in standard maneuvers to compare with those 
of a stable ship. A third objective was to apply DSL/360 
Language to simulate the ship motions. 
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II . EQUATION OF SHIP MOTION IN THE HORIZONTAL PLANE 



As developed in Ref. 2, the general motion of a ship is 
that of a rigid body in six degrees of freedom subject to 
gravity and buoyancy forces as well as controlling forces 
and hydrodynamic reactions and subject to hydrodynamic or 
other disturbances or excitations. 

Considering axes fixed to the body and parallel with 
the principal axes of inertia but with arbitrary origin as 
shown in Fig. 1, the linear and angular momentum equations 
may be written as 





( 1 ) 



where 



m = mass of ship 
U = linear velocity vector 



= ni + vj + w k 



u 



rate of surging 



v = rate of swaying 



r = rate of heaving 



= angular velocity vector 



= pt + qj + rk 



( 2 ) 



p = rate of roll = 0 



q = rate of yaw = 9 



r 



rate of pitch = 
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F = force vector acting on the body 

= XT + Y J + Zl (3) 

X = hydrodynamic force on body along x-axis 

Y = hydrodynamic force on body along y-axis 

Z = hydrodynamic force on body along z-axis 

hi = moment vector acting on body about a 

fixed point 

= kT + uf + n"£ (4) 

K = rolling moment about x-axis 

M = pitching moment about the y-axis 

N = yawing moment about the z-axis 

H = angular moment vector 

= I* p 4 + Iy q + I-y r k ( 5 ) 

I = mass moment of inertia about the axis 
z 

x, y and z respectively. 




FIGURE 1. Body axes 
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Inserting Equations (2), (3), (4) and (5) in Eq . (1) 

gives : 

X = 7 rj J u-t £ur-rv--*<r( f+r 1 ) T fc ( P%--r)+Zf(pr*pJ 
z = Z?\lirffir- $ ) + Zf(rp-f-)+^(rf-tp)} 

K = I x pl(Ii-fyp + m{faCw + fv--fu.)-%(v+ru.-pu;)} 
M =i2j+(Tj-T?)rp+}nl%(u--h$v- r --rr)-?i 4 (uri-f>v--l"-)} 

n "T^r-td^-l^Fl + Wj^fv + ru-pi^- fc(u. + fw-rv)} 



> (e) 




> ( 7 ) 



y 



In Eq . (6) and (7), the left hand side represents the 
forces and moments acting on the ship and the right hand 
side represents the dynamic responses of the ship, 
respectively. The forces and moments acting on the ship 
depend on the geometry of the hull, the motion and the fluid. 
That is, 

F (or hi ) = f (properties of ship, properties of 
fluid, properties of ship motion) 



Since the properties of ship and fluid are constant for 
a given ship in a given fluid without excitation forces, the 
forces or moments may be expressed as a function of the 
motion as, 



F(n. l?i) — j (%o y y 9j £ e/ j> / 0 / u. y ir, ur y p; r 

v, yp a h r, s , £, Jc. ) 



( 8 ) 
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where x q , y Q and z q represent a system of reference axes 
fixed in the earth and 






= orientation parameters 



o., 



v, ox, T, r I 

<r, or, P, $,r [ 

s, i, / 



motion parameters 



control surface parameters 



Considering motion in unrestricted calm water, only three 
degrees of freedom are of concern: yaw, sway and surge. 
Therefore, the forces and moments are independent of the 
change in orientation with respect to the earth axes. 

Furthermore since the forces and moments produced on the 

• . * 

ship as a result of £ and ^ are normally negligible. 

F( or /??) = f ( u , v, r, u, v, r , £) (9) 



Since the motion is restrained to the horizontal plane 
for the surface ship in calm water as prescribed in the 
previous section, roll, pitching and heaving are taken as 
zero, i.e., p=q=w=0 in Eqs . (6) and (7). Furthermore, most 

ships have = 0 so that Eqs. (6) and (7) may be written as, 



x = Tfifu -r\r~'X^T x ) 



y= i7i ( v - f- r u + Xq x ) 

N = I?r+^(^(v+ur)} 



( 10 ) 
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From Eqs . (9) and (10) then, these become: 

x: Tn(u-rv-Xtr') = JrKwr, v, r, S ) 

y: Pilvtra-v^r) =■ v, r, S) 

N . l,,r*M<,(v+ULr) - f„ ( a , V, K U,v, r, S) 



(ii) 
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Ill . 



DYNAMIC STABILITY OF THE SHIP 



Dynamic stability is directly related to the magnitude 
of yaw and sway deviations caused by small initial dis- 
turbances. In many problems related to dynamic stability 
with small disturbances, the linear analysis of motion is 
quite useful and will be applied in the following to the 
present problem. 

A. LINEAR EQUATION OF SHIP MOTION 

To simplify the equation of ship motion, the following 
assumptions are made: 

(1) The ship has adequate stability. 

(2) The water is so calm that the rolling, pitching 
and heaving motions are negligible. Therefore, the motion 
of ship may be regared as taking place in the horizontal 
plane only. 

(3) Yaw and sway do not affect the forward speed 
appreciably, and small changes in forward speed do not affect 
yaw and sway motions so that the surge equation is decoupled. 

(4) The ship has one rudder at the stern along the 
centerline and the effect of the propeller is disregarded. 

(5) The forces and moments on hull and rudder are 
expressed by means of velocity and acceleration derivatives. 

In order to linearize the equations of ship motion about 
an initial equilibrium • condition of motion with the previous 
assumptions, only the linear terms in the change of the value 
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of the variables from equilibrium, i.e. Au, Av, Ar, Au, av, 
Ar and are considered. 

From Eq . (9) then, 

X = X<>+ Xu At/ +X ir AV'-{ Xr A ? Xa i~Xir A V'+ A S 

Since in initial equilibrium V - To ~ u c ~ V — To — So - O . 

A\r^ IT- Vo = V 

A r = r-To = r 

au - a - u 0 - u 

A V - V - \To = v 

A r -r -To ~r 

^ S - S ~ So - s 

Therefore , 

X^AtX^lHXrV+X, r+/.u +XV+Xr r+Xi £ (12 > 

and likewise, 

y = % au + Yt^ y r+ /u i+ /^ + /f i+ YsS (l3) 

H + hiv 

where X , Y and N are the values of the initial equilibrium 
o o o 

conditions . 

In the initial equiliburium conditions of straight ahead 
motion at constant speed, there are no forces acting on the 
ship. Hence 

X = y. •= a/„ = ° 
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Furthermore, since the ship is symmetric about the x-z plane 



Yu = Y* =- X*= X r = Xf = hJu = Ncl = o 

and 

Xs^o 

Therefore , 

y* Xu AU + Xi a 

y=y v y r+Yr (r -t-y?X{-xrS o4> 

M = A /„ir+ A/, r + A/,/ V + A/r r+ A/r 

Linearizing Eq. (10) gives 

X -iii\ ( u 0 +AU)-(ro+*r)(i£i-Air)-%r( n-^r ) 2 } 

= M(AU-*rAr-KAr 2 ) 

This is simplified further by neglecting second order terms 
so that 

X- lr)Aa - mu 



and likewise, 

y s T/l( ¥+ JT u o + %f T ) 

% r ■+ (v + ru e ) 



(15) 



From Eqs. (14) and (15) 

( X ^- 711 ) 1 *+ Xu ALL 

( y* - /» ) IT + Kr V4 ( /f -W ; * ) ■ r + ( X - >» u. ) ■ r + X r = 0 

(^ <r -mXi)ir+hJ v - v '+(Xf-Ii)T+ Mr->uZfU,)r +^i=.o 



20 



For convenience, Eq . (16) is non-dimensionalized . Every 

term of the first and second equations for the X-force and 

the Y-force in Eq . (16) has the dimensions of force and the 

third equation for moments has the dimension of a moment. 

To non-dimensionalize these, the force equations are 

f a 

divided through by (-^-LHU) and the moment equation is divided 

p i 2 

by (T^-^^) where f= dinsity of ship, L = ship's length, 

U = ship speed and II = draft of ship. This gives, 





where 





for small perturbations, 
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In Eq. (17) the dimensionless time is defined as 



/'= -r* 



B. DYNAMIC STABILITY WITH FIXED CONTROL 

Taking the Laplace Transform of Eq. (17) with the rudder 
fixed at L- 0 gives: 

+ +(y r -to'uo)]rU)=o 

(iS) 

[o '?4)s+ \( r ]ic<;)A{(h/rIz)5 -HtJrft %W)J rU)=o 



Then, using Cramer's rule 



v(s) = 



_0 

A 



r'(s) 



0 

A 



where 



+ // 



A 






(/r -7ft %')s+(Y r - w) 

s / 

( /V/* ~Jjr ^(tJy ft) 'Xif ) 



Expanding this determinant gives the characteristic 
equation 



A=/1s-t-Bs-i*C 
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/ 



where 

A = (, Yv'~ w)(Hr -T?)-(, 'y/-^YY)(AJir-^% ) 

B =( /Y-MJfA/i -M 7 Y) ~/u- (A/r~Iz)- 

) ( 20 ) 

c = y v (A/f -W Xf )- h! \r ty* ~ to ) 



The roots of Eq. (19) are 



dT 



fo = 



-B -j /B-4/?C 






+ (&f-4& 



2 



( 21 ) 



B- _ fT& \ % 




-4& 



z 



Therefore, the free responses of ship, i.e., sway and 
yaw motions are given by 



/ art <r 2 t 

v(t')-'v,e -t-vie 

. <rf' fit' 

r ct') = r,e + r L e 



( 22 ) 



where v^, v 2 , r^ and r 2 are arbitrary constants depending on 
initial conditions. 

As shown in Eq . (22), the dynamical stability depends on 
the roots of the characteristic equation. Since (ff dz ^ 

JZf/t V ({') = 0 tZ/flcl JU*r> r'(t') -0 

t'-*CO / ’-i’CO 
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narginally 
uns taole 
rr.odel 8 



(i> 
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FIGURE 2. Response of the yaw rate and sway velocity to 



if and only if (J7 ^ 0 

/ / 

If then, yfoo)^ 0 and Y(O 0 )^O which means that the ship 

is dynamically unstable in straight line motion. The root, 
JT is called the 'stability index'. 

In Appendix B and computer program 1, a five foot Series 
60 model of the Davidson Laboratory with block coefficient, 
Cg = 0.7, and no propeller, was modified to be (a) stable 
(b) marginally unstable (c) unstable (d) very unstable and 
(e) extremely unstable. These models are referred to as 
model A, B, C, D and E, respectively. 

The hydrodynamic derivatives and stability indices of 
these models are in Table 1. 

Eqs. (21) and (22) were simulated to obtain the response 
of each model to small a disturbance. (See computer program 
7 and Fig. 2) 

As shown in Fig. 2 (which is the result of computer 
program 7), the yaw rate, r, as well as the sway velocity, 
v, of the stable model A decreases as time increases but 
the unstable models B, C, D and E tend to diverge due to 
the positive stability index. 



TABLE I. HYDRODYNAMIC DERIVATIVES AND 
STABILITY INDICES OF THE FIVE MODELS. 



iriod- 

& 


// 


a/; 


// 


V/ 


yY 




Yr' 


V/ 


(FT 


A 


-.?F4+ 




4.6774 


-.663$- 


-.mo 


+ . 004-0 


+ ,oo/y.o 


0/32 


-. 32/3 


B 


-•3/3/ 


"■//4F 


+-oir£7 




-'/B62 


-f .<3 05 //- 


■t.tfo ^ 4 - 


-.e/2-7 




C 


-.3/3/ 




+ .04-4-£ 


-.M-sr 




-f .0624- 


+.0 02.4- 


-.01 + 1 


+. 2^1 


D 


- ?/3' 




+ . oZ 6 5* 


-o3?o 


-. 126 % 


+.0010 


i-.oolo 


- .0117 


+. 01 % 


E 




- mi 


-.oUS 


— . OZI h 


— . tftoo 


0 . 


0 - 


-.oil?. 
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C. STABILIZATION OF THE DYNAMICALLY UNSTABLE SHIP WITH 
FIXED RUDDER 

As described in the previous section, a dynamically 
unstable ship has 




To stabilize this unstable ship, (fp must be made less 
than or equal to zero, i.e., (TT ^ 0 or, 




By the nature of hydrodynamic derivatives for the ship 
A^> 0 and B^> 0 always. Therefore, C must be greater than or 
equal to zero for stability. Thus, C is called the "stability 
criterion . " 

From Eq. (20) 

0 —Yxr (dr ~ W (/y -771 ) Y 0 (20a) 



Then, considering the characteristics of each of the 
hydrodynamic derivative, mass and the distance of the C.G. 
from the origin of the fixed coordinates on the ship. 

Y v is always large negative quantity 
N^ is always large negative quantity 
N v may be positive or negative 
Y^ is small but may be positive or negative 
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m is always large and positive 

Xq is always small but may be positive or 
negative 

Finally, N v has the greatest effect on dynamical stability 
since the term \ -W\) is a large negative quantity in 
equation (20a). If, for example, N v becomes more positive, 
the stability criterion C will become more positive so that 
the stability of the ship will be improved. It is noted 
that N v can be controlled by a lifting surface. 

1 . Sensitivity of Dynamic Stability to a Lifting 
Surface 

Equation (21) was simulated using the computer 
program DSL/360 Language for the effect of the location, size 
and aspect ratio of a lifting surface (movable fin) on the 
models A, B, C, D and E. (See the computer program 2). The 
results are summarized in Table II and III. 

The following are indicated by Tables II and III and by 
Figs. 3, 4 and 5: 

1) The position of the fin at which the stability index 
is maximum is independent of the area and aspect ratio of the 
fin, because this position is the pivot point at which the 
relative flow is parallel to the centerline of the ship. 

(See Fig. 14) 

The parallel flow makes no attack angle on the fin so 

that Y , N , Y and N of the fin are zero but Y‘ , N' , Y* and 
v v r r v v r 

NL ai'e not zero as described in Appendix A. Therefore, the 
stability criterion C in Eq . (20) is the same as that of the 
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TABLE II. EFFECT OF FIN AREA AND FIN LOCATION WITH 
FIXED ASPECT RATIO ON DYNAMICAL STABILITY 



Fin 


X f/L 


A 


B 


C 


D 


E 


No fin 




-0.3213 


+0.0238 


+0.2561 


+0.5194 


+1 . 3438 


0 . ..=1.5 
ef f 


- 0.5 


- .4366 


- .1066 


+ .1108 


+ .3625 


+1 . 1684 


A f =0.0056 


m 


- .3492 


- .0020 


- .2324 


+ .4968 


+1 . 3223 


h f =0. 015 


pivot 


- .32079 


+ .0242 


'+ .2569 


+ .5197 


+1.3441 


point 


(0.48) 


(0.4) 


(0.34) 


(1.305) 


(0 . 27) 




+ 0.5 


- .3208 


+ .0226 


+ .2512 


+ .5109 


+1.3288 


d =1.5 
ef f 


- 0.5 


- .5707 


- .2616 


- .0642 


+ .1710 


+ .9474 


A f =. 01265 


a 


- .3834 


- .0341 


+ .2026 


+ .4682 


+1.2944 


h f =.0225 


pivot 


- .3203 


+ .0246 


+ .2567 


+ .5199 


+1.3442 


point 


(0.48) 


(0.4) 


(0.34) 


(0 . 305) 


(0.27) 




+ 0.5 


- .3204 


+ .0212 


+ .2458 


+ .5011 


+1.3110 




- 0.5 


- .7376 


- .4583 


- .2887 


- .0780 


- .6717 


A f =.0225 


£ 


- .4299 


- .0780 


+ .1617 


+ .4286 


+1.2548 


h f =.03 


pivot 


- .3197 


+ .0250 


+ .2571 


+ .5202 


+1 . 3444 


point 


. (0.48) 


(0.4) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


- .3199 


+ .0196 


+ .2393 


+ .4893 


+1.2888 


d --=1.5 
ef f 


- 0.5 


-1.0956 


- .8930 


- .7921 


- .6486 


- .0573 


A f =0.0506 


& 


- .5532 


- .1963 


+ .0507 


- .3198 


+1 . 1427 


h f =0.045 


pivot 


- .3185 


+ .0260 


+ .2579 


+ .5209 


+1.3448 


point 


(0.48) 


(0.4) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


- .3188 


+ .0162 


+ .2259 


+ .4643 


+1.2393 


destabilizing 

area 


0.415i-~ 

0.5L 


0. 35L~ 
0.445L 


0. 3L — 
0. 375L 


0 . 27L 
0.34L 


0.24L- 
0. 3L 



Note: All values in column A,B,C,D and E (except the values 

in ( ) which mean pivot point distance fro-m $£) are 

(J^ of model A , B , C , D and E, respectively. 
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TABLE III. EFFECT OF FIN LOCATION AND ASPECT RATIO WITH 
FIXED AREA ON DYNAMICAL STABILITY 



Fin Type 


X f/L 


Model A 


B 


C 


D 


E 


No fin 




-0.3213 


+0.0238 


+0.2561 


+0.5194 


+1.3438 


%/= 1.5 


- 0.5 


-0.4366 


-0.1066 


+0. 1 108 


+0.3625 


+1.1684 


A f =0.0056 


£ 


-0.3492 


-0.0020 


+0 . 2324 


+0.4968 


+1.3223 


h f =0. 01 5 


pivot 


-0.3208 


+0.0242 


+0.2569 


+0.5197 


+1.3441 


point 


(0.48) 


(0.40) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


-0.3208 


+0.0226 


+0.2512 


+0.5109 


+1.3288 


&jf= 2. 


- 0.5 


-0.4557 


-0.1275 


+0.0872 


+0.3373 


+1 . 1442 


A f =0.0056 


& 


-0.3539 


-0.0063 


+0.2285 


+0.4932 


+1.3189 


h f =0. 0173 


pivot 


-0.3207 


+0.0242 


+0.2565 


+0.5197 


+1.3439 


point 


(0.48) 


(0.4) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


+0.3208 


+0.0224 


+0.2505 


+0.0510 


+1.3265 


% 3. 


- 0.5 


-0.4823 


-0.1572 


+0.0543 


+0.3021 


+1 . 1036 


A f =0. 0056 




-0.3605 


-0.0122 


+0.2231 


+0.4881 


+1 . 3164 


h f =0.0212 


pivot 


-0.3206 


+0.0242 


+0.2565 


+0.5197 


+1 . 3440 


point 


(0.48) 


(0.40) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


-0.3207 


+0.0221 


+0.2495 


+0.5078 


+1 . 3235 


II 


- 0.5 


-0.5000 


-0.1769 


+0.0325 


+0.2788 


+1.0786 


A f =0.0056 


$ 


-0.3649 


-0.0162 


+0.2195 


+0.4847 


+1.3112 


h f =0. 02448 


pivot 


-0.3206 


+0.0244 


+0.2565 


+0.5198 


+1.3441 


point 


(0.48) 


(0.4) 


(0.34) 


(0.305) 


(0.27) 




+ 0.5 


-0.3206 


+0.0219 


+0 . 2488 


+0.5066 


+1.3215 



Note: All values in column A,B,C,D and E (except the values 

in ( ) which mean pivot point distance from J£) are 

of model A,B,C,D and E, respectively. 
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fin position 

3 The effect of fin location and fin area with 
constant aspect ratio on stability of the ship. 
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Fig. 




constant area of fin on the stability of the ship. 
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upper limit 



A 
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Fig. 5 Destabilizing area of the fin and the most destabilizing point 



no-fin condition. In other words, the stability criterion, 

C is minimum at the pivot point during shifting of the fin 
from stern to bow. 

This pivot point is the position of the fin at which it 
exerts maximum destabilizing influences on the ship. 

2) The effect of the fin on stability increases with 
fin size. 

3) The more unstable the ship, the closer the pivot 
point moves toward the midship. 

4) Compared to the no-fin condition, the placement of 
the fin near the pivot point has a negative effect on 
stability. This criterion is called the "destabilizing area" 
in this thesis. This destabilizing area is due to the 
relative change in A, B, and C in Eq. (21) as the fin location 
is shifted. C is almost constant near the pivot point and 
the destabilizing area j.s independent of the size and type 

of the fin. 

5) Moving the fin toward the stern from midship 
improves the dynamic stability of the unstable ship. Con- 
versely, moving the fin forward of midship does not cause 
much change from no-fin condition regai’dless of the size or 
type of the fin. 

6) Stability is sensitive to the fin area but rather 
insensitive to changes in the aspect ratio of the fin. 

2 . Sensitivity of Stability of the Unstable Ship to the 
Position of the C.G. of the Ship 

It is clear that moving the C.G. toward the bow tends 
to make 'C' more positive in Eq . (20). The effect of moving 
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C.G. was studied by use of Eqs . (20) and (21) using DSL/360 
Language for model C. The results are shown in Table IV. 
(See computer program 3). 

i 

TABLE IV. THE EFFECT OF LOCATION OF C.G. ON 
THE STABILITY INDEX 



X G/L 


-0.4 


-0.25 


& 


+0.05 


+0.1 


+ 0.25 


6 


+ .6424 


+ .4022 


+ .0775 


+ .015 


- .0473 


- .2297 



As shown in the above table, stability was improved by 
moving the C.G. toward the bow as expected. However, this 
is not a practical method of stabilizing an unstable ship. 

D. STABILIZATION OF THE DYNAMICALLY UNSTABLE SHIP WITH 
RUDDER CONTROL 

In previous sections, the stabilization of a dynamically 
unstable ship using lifting surface was discussed. However, 
it is clear that the dynamically unstable ship can be made 
directionally stable if the rudder can be forced to compen- 
sate for disturbances so that the ship will maintain course 
with only small oscillations in yaw and sway. If the 
deflection of the rudder,^, is proportional to the heading 
error, ~\jf , and the angular velocity, r = 'j r > then the equation 
of rudder deflection becomes 

&({) = k ir(t') + kii'tt’) (23) 
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where and/?^ are constants of proportionality of the 
control system. 

Inserting Eq. (23) into Eq . (17) and taking Laplace 
Transform gives 




( y r -m'+Y/kp ip k] f^>o 



(24) 




As these algebraic equations of (24) are homogeneous, a 
non-trivial solution (v \ *'/' ^ 0) demands the determinant of 
the system, i.e., characteristic equation, to be zero. 

Let the characteristic equation be 

A 5 3 + B s* + C s + V = o ( 25) 



where 



3 =(y^-m')(A/r^y^-y ds) + X-V/-I/; - 

A /„' ( 'Yi- m '>Y) - ( y r '-»t '+ k /{, )(Yir ~ m ^ ) 

c = )- 

blU y m '+ 4 y/ ) - 4 //( a 4 -> n V ) 



/ , / 



d ~ k, yy^s - k % 
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(T 3 • Then, 

(26) 

* f f I 

If and only if dT , (T\ and <T$ are negative, then y"( oa) and 

f 

V(oo) will be zero and the unstable ship will be stable. 

Furthermore, if the roots are three repeated roots, i.e., 

/ / / 

(5T = (Tz = (Tj <( 0, the ship will be stable with the shortest 
possible response. 

By designing for the optimum value of and in the 
control mechanism, the ship will be automatically controlled 
for suitable stability. 

To find the three repeated negative roots of the 
characteristic Eq. (25), the 'MITROVIC METHOD' was applied. 
(See Appendix C). 

In computer program 4, the optimum k, and for the three 
repeated negative roots of Eq . (25) were calculated for model 

C and Eq. (26) was simulated by the DSL/360 Language. (See 
computer program 4) . 

As shown in Fig. 5-1 which is the result of computer 
program 4, the response with proper k, and k z for three 
negative real roots of characteristic equation (25) is an 
exponential transient. Conversely, if the auto-steering 
device has some k t and f? z which give complex roots of charac- 
teristic equation, the response will be stable with some 
oscillations . 



/ / 



Let the roots of Eq . (25) be (dT dz anc * 






/ 07 k' Gi't' Jl't' 

V(i')- v; e + 
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The auto-steering device should be able to select the 
proper control costants, k t and /?■>_ . 

It is evident that the dynamically unstable ship can 
be made to be stable directionally with auto-steering control. 
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_L 

lO 




(u'btp'c^t) Suxpueq s.diqs 
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Figure 5-1. Response of ship's heading by auto-steering control after small disturbance. 



IV. MANEUVERABILITY OF THE UNSTABLE SHIP 



To predict the performance of a dynamically unstable ship 
as well as to get realistic predictions of maneuvers with 
large rudder angles, a non linear mathematical model is 
needed . 

A. NON-LINEAR EQUATION OF SHIP MOTION. 

In the non-linear analysis of the ship motion, the 
following assumptions which are rather typical in the 
literatures Refs. 2 and 8 were applied. 

1) Terms of higher order than third order are neglected. 

2) Rudder force and moment derivatives of higher than 
first order, and effects of rudder angular rate are negligible. 

3) Odd terms in v, r, & , v and r in the X-force equation 
and even terms in the same parameters of the Y-force and 

N-moment equations are eliminated due to symmetry. 

4) Cross coupling terms between acceleration and 
velocity are neglected. Taking the Taylor Series expansion 
of Eq. (11) with the above assumptions: 



. / 




U - 
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r = 



'-//) i(u' v, r', Y)- (in % - U^) ( < v, r' Y ) 



where 



7 2 



( b'-y^d^r ') - (7n%'-K /y)(7r>'%'-yr ') 

i,(u,v, y'/ <0 =Xo+ Xu Au + 2 Xmu al( +t x “«u 4m + 

+ X Xf ) r +iA& 

^ V+ j- x ; ruT '^u' 4 

(Xvr + 7,1 7 l/r " X Xv^ V" X Xff ' Y J- XvnYy'^d 
Xv^u^ AR +^r^u r ^ AC(/ z 

£.(u,v / r', Y)- X+Xu 4U +X uuAU + Xr r + X’X^ v "x 

)^ tt XkXX/juu^ 11 -i )/ r ^Vv' / X 
f ? ( u'v;' r' <$) = a/„V A/o« + M> uu + a/ v v + s'A/ mr v + 

~^bivxr^ ^ ^ ^ ^"z^irUU 

6^r “ w V+ r rr ^ + i: a//^ 1 v ' ^jrri'ru r % *+ 



A / ru Y*u'+z^ruu Y*U' 4- A/^ i + ^ 

T^iro -!■■£■ fJ fa &’ + hl£ u £ Aa "+ 

-^dJ UIJ Y^x dy ri vrd 



(28) 
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Integrating Eq . (27) 




( 29 ) 



o 



Note: The time domain should be non-dimensional. 

So far the reference axes x, y and z are fixed to the 

moving ship. But it is convenient to take the reference axes 

x and y fixed to the earth in order to determine the path 
o o 

of the ship. Adopting the figure from Ref. 1 gives 




— x, 



FIGURE 6. Relations between axes fixed 



in the earth and axes fixed in the ship. 
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The following relationships may be written as 



X(i) = u (i) a* - v(t) /tin V ' (t) 



( 30 ) 



S O 



( 31 ) 



Itf) = f r i (t) dt 



Non-dimensionalizing Eqs . (30) and (31) using the 

del' initions. 






u 



X 



u 









it 



u 



/ 



/ u 
t =/ '[_ 





( 32 ) 



Likewise , 
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B. TURNING CIRCLE WITH SELECTED RUDDER DEFLECTION 



Fig. 7 shows the typical rudder deflection with respect 



to time in a turning maneuver. 




FIGURE 7. Rudder deflection vs. 
time in turning maneuver 






Using the ramp function for rudder deflection, then 



j ({) = o' [ ranfCk-tuo) - 



(33) 



where 

S> = rudder deflection 
ramp (t) = ramp function 

^ = rudder deflection rate 

= time delay 

^<W = time at Smart , 

Calculation of steady turning radius. 

The relation of the angular velocity with linear 
velocity gives : 

(radius) x (angular velocity) = (linear velocity) 
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Ult) J~U Z (I)+U ? (t) 

or Radius = — 

mi m> 



Taking non-dimensional terms 



’^o.cSJuUl. 



L 



U /l/V V x 

U L r 






R&d 




(34) 



where 



JdaA 



RaJ-dtuL 

_ 



Figure 8 shows the typical turning trajectory and 
terminologies corresponding to Table 5. Eqs . (31) and (34) 

were simulated for the 5 models using computer program 5 in 
turning circle maneuvers with a rudder deflection rate of 
2 degrees per second for the prototype (L = 500') with zero 
time lag and constant final angles of 10°, 20° and 35°. 

The results are presented in Table 5 and Figs. 8a, 8b and 8c. 



FIGURE 8. 
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-1.00 O.TO 2.00 2. TO S. DO <5. TO 







-2.5 



“KT 

I 



.5 



1 . 



2.5 



4. 



5.5 



7. 



8.5 



X / L 
o ' 



FIGURE 8a. Turning circle of five models with 
35 degrees of rudder angle. 
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curve # 1 : model A> #2 : model B, #3: mo del 3, model D 
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\n 
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5 1. 2.54. 5.5 7. 8. 

X Q / L 

FIGURE 8b. Turning circle of five models with 
25 degrees of rudder angle. 
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curve /^lsirodel A 3 //2:ro lel h, ^3:noiel 3, #U:roiel D, #£;rroisl E 




! | | | | — | 
1.6 3.6 5.6 7.6 9.6 11.6 

X c / L 

FIGURE 8c. Turning circle of five models with 
10 degrees of rudder angle. 



47 






TABLE V. CHARACTERISTICS OF MODEL A,B,C,D AND E 

IN TURNING CIRCLE 





Model 


Adv. 


Max . 
Adv. 


Tran . 


T.D. 


Max . 
Tran . 


S.T.R. 


*'/ 

' rev . 




A 


6.86 


6.8825 


4. 86 


10.30 


10.322 


5.0722 


31.6 




B 


5.30 


5.3224 


3.25 


7.08 


7. 1079 


3.4498 


21.4 


10° 


C 


4.49 


4.5088 


2.65 


5.88 


5.7079 


2.7611 


17.2 




D 


3.94 


3.9667 


2.18 


4.79 


4.8122 


2.3290 


14 . 6 




E 


3.08 


3.1277 


1.45 


3.51 


3.5490 


1.7295 


10.8 




A 


4.80 


4 . 8269 


3.16 


6.91 


6.9411 


3.4177 


21.2 




B 


4.07 


4.0996 


2.40 


5.46 


5.4962 


2.6946 


16.8 


20° 


C 


3. 60 


3.6306 


2.06 


4.66 


4.6863 


2.2951 


14.2 




D 


3.29 


3.3134 


1.80 


4.10 


4.1379 


2.0266 


12.6 




E 


2.70 


2.7892 


1.35 


3.23 


3.2719 


1.6070 


10.0 




A 


3.840 


3.8685 


2.33 


5.2203 


5.2433 


2.6364 


16.2 




B 


3.430 


3.4584 


1.98 


4.435 


4.4551 


2.1938 


13.6 


35° 


C 


3. 120 


3.1538 


1.77 


3.900 


3.9380 


1.9385 


12.0 




D 


2.918 


2.9475 


1.49 


3.540 


3.5804 


1.760 


10.8 




E 


2.550 


2.5943 


1.15 


2.790 


2.9890 


1.4727 


8.6 



Note: Adv. = Advance/L 

Tran. = Transfer/L 

T.D. = Tactical Diameter/L 

S.T.R. = Steady turn radius/L 
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Fig. 9 Relation's between turnirg circle paraneters and 

stability index at 3*? degree of rudder deflection 
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Fig. 10 Advance vs. rudder angle in turning circle. 
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Fig# 11 Tactical diameter vs. Rudder angle in turning circle# 
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^ig» 12 Steady turning radius vs. rudder angle 
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'ig.13 Circling tire vs. rudder argle in turning circle. 



Discussion of turning circle maneuver 

1) As shown in Fig. 9, the dynamically less stable ship 
has a smaller tactical diameter, advance, steady turning 
radius and circling time than the more stable ship, because 
the yawing rate of the unstable ship is greater than that of 
stable ship at the same rudder deflection. However, the 
turning circle parameters are becoming constant as the 
stability index increases. 

2) As shown in Figs. 10, 11, 12 and 13, the rate of 
change of the turning circle parameters versus the rudder 
deflection is high at the small rudder angle. In other 
words, the above parameters are less sensitive to rudder 
angle as the rudder deflection increases. 

3) The more unstable the ship, the less sensitive the 
above parameters are to the rudder deflection. 

1 . The Effect of the Location of the Rudder on the 

Turning Circle 

A change in the location of the rudder has multiple 
effects on the ship motion. One is it's effect on dynamical 
stability, which is the same as the effect of a movable fin 
described in section (II-C-1) and the other is it's effect 
on maneuverability. 

As described in Appendix A and Ref. 1, 

1) Y , and Y<$ are independent of the position of 
the rudder. 

2) N , Y , N* and-Y- are less positive if the rudder 
moves forward from the stern. 
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3) N r and are less negative if the rudder moves from 
the stern to the midship position and less positive from 
midship to the bow. 

4) Since the relative flow passing the ship depends on 
the position considered, the attack angle depends on the 
location of the rudder. Forward of the pivot point the 
attack angle is larger than rudder angle, at the pivot 
point the attack angle is equal to rudder angle and aft of 
the pivot point the attack angle is smaller than the rudder 
angle . 

5) N c c strongly depends on the location of the rudder 

since -Xf where is the distance of the rudder 

position from midship and is negative aft of midship. N c ' 
is less negative if the rudder moves forward to the bow, 
but N^o is less positive if c is negative. This less 
positive effect decreases the yawing rate, r, due to 
decreases of f^ in Eqs . (23) and (28). Decreasing r in 

turn causes the turning radius to increases. 

On the other hand, the movable fin effect of shifting 
the rudder forward makes the ship more unstable which makes 
the turning radius decrease. These two compromising effects 
are shown in Fig. 14 and 15 which were obtained from 
computer program 6. In computer program 6, model C was 
simulated for the evaluation of the effect of the rudder 
position on turning maneuvers. 

As shown on Fig. 14 which represents the effect of the 
rudder position on the turning radius without replacement of 
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TABLE VI. THE EFFECT OF RUDDER POSITION ON THE TURNING 
CIRCLE OF MODEL C WITHOUT A FIN AT THE STERN. 





t' // uur . 


v' 


/ 

r 


RaJL ' 


€ 


- 0.5 


7 . 12 


- 0.1617 


0.4413 


2.2960 


0.25616 


- 0.4 


6.84 


- 0.16571 


0.45764 


2.25 




- 0 . 3 


6.72 


- 0.1701 


0.4688 


2. 164 


0.5482 


- 0.22 


6.64 


- 0.1735 


0.4734 


2.144 


0.6492 


- 0.2 


6.6 


- 0.1743 


0.4739 


2 . 14 




- 0.17 


6.64 


- 0.17556 


0.4742 


2.1412 


0.70581 


- 0.10 


6.64 


- 0.1784 


0.424 


2.15 




% 


6.76 


- 0 . 1822 


0.4639 


2.19 


0.8552 


+ 0.1 


7.0 


- 0.1856 


0.4485 


2.268 


0 . 90969 


0.15 


7.16 


- 0.1871 


0.4383 


2.32 




0.20 


7.44 


- 0.1885 


0.4265 


2.386 


0.9394 


0.21 


7.43 


- 0 . 18869 


0.42398 


2.4 


0.9410 


0.22 


7.45 


- 0 . 1889 


0.4214 


2.415 


0.9424 


0.24 


6.08 


+ 0.1345 


- 0.51396 


1.963 


0.9445 


0.25 


6.08 


+ 0 . 1352 


- 0.5163 


1.955 




0.27 


6.04 


+ 0.1366 


- 0 . 5205 


1.94 




0.30 


6.0 


+ 0 . 1385 


- 0. 5260 


1.92 


0.9456 


0.34 


5.92 


0. 14099 


- 0.53204 


1 . 898 


0.9420 


0.5 


5.8 


0.1493 


- 0.5425 


1 . 804 


0.8980 
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steady turning radius/L 




?ig.lh The effect of the position of the rudier or 
turning radius without a fin at the stern. 
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TABLE Via. THE EFFECT OF RUDDER POSITION ON 
THE TURNING CIRCLE WITH A STABILIZING FIN AT THE STERN. 



#)/l 


v' 


/ 

r 




d 7 


- 0.5 


- 0.1193 


0.3053 


3.299 


- 0.3804 


- 0.4 


- 0.1233 


0.3112 


3.238 


- 0.2570 


- 0.35 


- 0.1252 


0.3125 


3.225 


- 0.1937 


- 0.33 


- 0.1259 


0.3127 


3.223 


- 0.1683 


- 0.3 


- 0.1269 


0.3126 


3.225 


- 0.1305 


- 0.25 


- 0.1286 


0.3113 


3.238 


- 0.0625 


- 0.2 


- 0.13 


0.3086 


3.268 


- 0.0088 


- 0.1 


- 0.1322 


0.2983 


3.382 


+ 0.1001 


0 . 


- 0. 1333 


0.2810 


3 . 590 


+ 0.1892 


0.2 


- 0.1292 


0.2220 


4.542 


+ 0.2948 


0.24 


- 0.1267 


0.2054 


4.908 


+ 0 . 3039 


0.26 


- 0.1251 


0. 1961 


5 . 138 


+ 0.3070 


0.27 


- 0.1242 


0.1913 


5 . 268 


+ 0.308 


0.28 


+ 0.0572 


- 0.2827 


3.544 


+ 0.309 


0.30 


+ 0.0621 


- 0.2971 


3.373 


+ 0.310 


0.4 


+ 0.0779 


- 0.3413 


2.939 


+ 0.305 


0.5 


+ 0.0874 


- 0.3639 


2.758 


+ 0.2820 




- 0.1617 


0.4413 


2.296 


+ 0.2561 



Turning radius 




Rudder position / L 

FIGURE 15. The effect of the rudder position on turning 
circle of model C with a stabilizing fin at the stern. 
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Stability index 



the fin at the original rudder location, the steady turning 
radius decreased as the rudder was moved from -0 . 5L to 
-0.2L, and increased as the rudder was moved from -0 . 2L to 
+0.22L. At the rudder position 0.24L, the turning radius 
decreased suddenly from 2.415L to 1.963L. Also, since the 
turning moment at this point was reversed in direction, the 
ship suddenly swings to the opposite side. This phenomenon 
is the characteristic of the unstable ship as described in 
section IV-D. 

If the turning moment is decreased due to shifting the 

1 

position of the rudder, the yawing rate decreases until the 
rudder change occurs. After changing the turning direction, 
the magnitude of yaw rate was suddenly increased. Finally, 
the turning radius was decreased markedly. 

To find the optimum position of the rudder for proper 
stability and maneuverability of model C, the rudder was 
moved forward to the bow and in it's location a stabilizing 
fin was installed. 

As shown in Fig. 15 and Table Via, the turning radius 
at rudder position -0 . 33L is a minimum before a sudden 
change of rotation and in addition the stability index is 
negative. Although this minimum turning radius (3.223L) is 
greater than that of the original model C (2.296L), it is 
an acceptable optimum rudder position due to it's better 
stability index, i.e., negative index. 
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C. ZIG-ZAG MANEUVER 



According to the standard zig-zag maneuver, the relation 
of the rudder deflection ( £ ) and ship's heading ( *'/'“) is as 
described in Fig. 16. 




ship's heading and time in standard zig-zag maneuver. 



Using a ramp function for the rudder displacement, the 
rudder angle is given by : 

Sd) = cx[mrf {f'-i, ) - mpit-h ) - rmp(t-h ) + 

Y o/ro,f (x -t v )r Tmf ( t- Xg ) ] 
for / cyck utiXtf A< 4 , 

where 



u = 


rudder 


rate ( degrec/sec. ) 


Ympit)= 


ramp 


function 


t, = 


time 


delay of rudder 


/; = 


time 


at 


6- Cn-AX 


A = 


time 


at 


cf ='f 


ii 


time 


at 


<Y- ~ c rrax 


ii 


time 


at 


l=t 
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FIGURE 18. Yawing Rate and Sway Velocity in Z-maneuver. 
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FIGURE 19. 



Path Trajectory in Z-maneuver. 
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In computer program 8, the standard zig-zag maneuvers 
were simulated for model A and C. As shown in Figs. 18, 19 
and 20 which represent the results of the program 8, the 
following can be observed: 

1) The amplitude of overshoot in yaw angle of the 
unstable ship is larger than that of the stable ship. 

2) Although the first reversing execution of the rudder 
of the unstable ship is earlier than for the stable ship, 
the period of the stable ship is shorter than that of the 
unstable ship. 

3) Amplitude of overshoot of path in the unstable ship 
is greater than stable ship. 

4) The reason for 1), 2) and 3) is that the yaw rate of 
the unstalbe ship is greater than that of the stable ship but 
the sway velocities are almost the same for either ship so 
that it takes longer for the unstable ship to get the reverse 
heading after reversing the rudder angle. 

D. SPIRAL MANEUVER 

The standard spiral maneuver of model C was carried out 
using DSL/360 in computer program 9. The process was 
started with 10 degrees port rudder and 1 degree of deflect- 
ing increment and ended with 10 degrees starboard rudder. 

The process was repeated. 

Figure 20 and Table 7 show the results of computer 
program 9. 

During the spiral test of model C, the angular velocity 
was reduced suddenly from r' = 0.19 to r' = 0.28 at £ = 1.5 
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FIGURE 20. Spiral Maneuver of Model C 
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TABLE VII. SPIRAL MANEUVER OF MODEL C 



RUDDER ANGLE 
(DEGREE) 

******* 

- 9.000 

- 8.000 

- 7.000 

- 6.000 

- 5.000 

- 4.000 

- 3.000 

- 2.000 
- 1.000 
- 0.000 

1.000 

2.000 

3.000 

4.000 

5.000 

6.000 

7.000 

8.000 

9.000 

10.000 

11.000 

12.000 

11.000 

10.000 

9.000 

8.000 

7.000 

6.000 

5.000 
• 4.000 

3.000 

2.000 

1.000 

- 0.000 
- 1.000 
- 2.000 

- 3.000 

- 4.000 

- 5.000 

- 6.000 

- 7.000 

- 8.000 
- 9.000 



YAW RATE 
(RAD. /SEC*L/U) 

0.36548072 
0.35614842 
0.34629613 
0.33535161 
0.32470226 
0.31269461 
0.29961270 
0.28513420 
O’. 26874071 
0.24949914 
0.22536594 
0 . 18941998 
-0.28294164 
-0.29785007 
-0.31125194 
-0.32350868 
-0.33485866 
-0.34546858 
-0.35545886 
-0.36492395 
-0.37393415 
-0.38254720 
-0. 39080832 
-0.38356121 
-0.37394869 
-0.36493850 
-0.35547382 
-0.34548380 
-0.33437481 
-0.32352638 
-0.31127167 
-0.29787230 
-0.28296709 
-0.26596129 
-0.24573880 
-0.21966898 
-0.17531638 
0 . 28510946 
0.29959112 
0.31267542 
0.32468486 
0.33583581 
0.34628141 
0.35613346 
0.36548072 
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degrees port side. That means any increase of the rudder angle 
to port by 1.5 degrees will cause the model C to suddenly 
swing from starboard to port until r' = 0.28. Likewise, when 
the rudder angle increase from port to starboard the model C 
can turn against its rudder byrf= 2.5 degrees starboard and 
then suddenly swing in the opposite direction to a new stable 
position r' = 0.28. 
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V. CONCLUSION 



This thesis has presented the characteristics of the 
dynamically unstable ship in stability and maneuverability . 

As discussed in Part III, the most powerful factor in 
stabilizing a unstable ship is a lifting surface installed aft 
of midship; the addition of area of fin is more effective 
for stability than the aspect ratio. The effect of a lifting 
surface may be obtained by designing a fine stern with neat 
flow lines and deadwood or a proper rudder. It was found 
that the fin installed near the pivot point destabilizes 
the ship. 

As discussed in part IV, the unstable ship is more 
maneuverable in a turning maneuver but less maneuverable in 
the standard zig-zag maneuver than the stable ship. In 
turning maneuver, the more unstable ship was less sensitive 
to the rudder deflection. 

As discussed in section IV-B-1, the followings are 
considered necessary for optimum stability and maneuverability 
of a ship. 

1) Locating the rudder at the point aft of the midship 
where the steady turning radius is minimum. 

2) Installing a stabilizing fin or deadwood as the 
stern in order to maintain adequate stability. 

DSL/360 Language was found to be very powerful for 
simulating ship motion. 



69 



RECOMMENDATION 



VI . 

All results in this thesis were obtained by simulation 
using DSL/360 without real model tests. 

The author would like to compare these results from 
real model tests in the future. The results in the thesis 
are in non-dimensional values so tha-t they can be converted 
to full scale ships or models. 

Also the following itmes are recommended for future study; 

1) Maneuverability of the unstable ship in Z-maneuvers 
other than the standard Z-maneuver. 

2) Maneuverability of a dynamically unstable ship with 
auto-steering control. 

3) Ship design for the optimum stability and maneuver- 
ability . 
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APPENDIX A 



Prediction of Hydrodynamic Coefficients in Linearized Equa- 
tion of the Ship. 



Using the method in Ref. 1 and 2, the following predic- 
tions can be obtained. 

(1) A fin travelling at a forward velocity u and trans- 
verse velocity v has 




where 



A f / L H 

Aj rr ftyjlCL ojf 

(3 = 

9 Cl 



2 71 



- 711 a PjJ 



2 



H 
u it 



for general aspect ratio fins 



for very low aspect ratio fins 




(34a) 



-CXf'fWf 
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/ 



(35) 






f 



2irti t A{ 
(a + \)' A 



where 



a. 



aspect ratio 




effective aspect ratio 

Vc 

effective span 




(35a) 



Considering the rudder as a fin, 

y, <—K>{ 



(35b) 
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(2) Considering the ship's hull as a fin with consider- 



able drag coefficient C 




(36) 



where 




7T n - 



IT H 
L 



(l\!v ~ ' k, w ') * ( /v 

( Yx )/i - -A ^ ■+ ( yj )h 



( ^r)h *** ^ /"v )h 

Ov ^ “ ~^ ? 2. 




(36a) 



^ A// ^ ~ 0 



(see Ref. 1 for more details) 
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(3) Assuming a large deadwood which has a sufficiently 
low aspect ratio, an effective span equal to the draft of the 
ship and a triangle shape, than 




(37) 
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(4) Assuming that a rudder or small deadwood has 



'ZCu 



2.T 




due to the high aspect ratio, then 




(see Ref. 1 and 2 for more details) 
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APPENDIX B 



Models for Stable, Marginally Unstable, Unstable, very 
Unstable and Extremely Unstable Ship . 

In this section it is intended to obtain models for 
several types of unstable ships. These will be simulated 
using DSL/360 to study their stability and maneuverability. 
The original model was adopted from a five foot Series 60 
model (with block coefficient 0.7 and no propeller) as 
originally developed by the Davidson Laboratory. 



Principal characteristics of parent model. 



Length (L) 

Breadth (B) 

Draft (H) 

Block coefficient (C ) 

LCG from bow (XA 

(j 

Displacement (A) 

Area of rudder (A n ) 

Rudder span 

Rudder chord 

Mass (m) 

Mass coef f . (m') 

Longitudinal added mass 

coef f . (k-^m ' ) 

Lateral added mass coeff. (m' ) 

2 



5.0 ft. 

0.714 ft. 
0.267 ft. 

0.7 

5.15 ft. 

41.64 lbs. 
0.021 sq . ft. 
0.2 ft. 

0.105 ft. 
1.292 slugs 
0.2 

0.004 

0.18 
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Rotational added mass coeff. (m') 0.165 

2 

Radius of gyration in air 0.25L 

C. G. of lateral added mass from 

midship (— ) 0.024 

Drag coeff. at zero drift angle 0.019 



By using Appendix A and computer program 1 , the parent 
model can be modified to the desired 'models . From the 
results of computer program 1, the five models were selected, 
which were named stable model A, marginally unstable model B, 
unstable model C, very unstable model D and extremely un- 
stable model E. (See computer program 1, Appendix A and 
Table I) 

The hydrodynamic derivatives in Table I were considered 
to be inherent for each model in this thesis. 
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APPENDIX C 



Selection of the Optimum Values of the Proportionalities 
of the Rudder Deflection to the Ship's Heading error and 
Yawing Velocity . 

Let the rudder deflection equation be 

Sd) = /,'>/' + <see E <i- 23) 



The ship motion under such rudder control has the 
characteristic equation, 



A s -v 3 s C $ P ^ 0 



(see Eq . 25) 



Let 



B „ 

■s = -a r 



, then 



B 



T F’ •+ — — - P -f- p -f V - O 



C 3 



Pi 



A 



P 



or 



5 2 AC A 2 L> 

P + P + 4 s P + e* : =0 



or 



P 3 A- p 4 P + b 0 - 0 



(39) 



) 
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where 




and 




(40) 



From Eqs . (24) and (25) 




3 — 2 + C2. 



(41) 



wher’e 







Likewise 



where 



C = § A + c 4 + c r 

% = ^s ( Vi 

c- %') - /vyiy/V 



(42) 



V ' V/ 



(43) 



where 



,/ , / / 



^ _ ( A/cj A/y x? 

Eq. (41) gives 

b % =c;^+2c,c/^c 

S 5 = C ; -h 2 C^c,. hi + 2 c, g 4 + C, 3 

Eq. (40) gives 

A£_ 
b> ;= 

4 ^ 2 4V^C/ c 2 4+ c 4) M (' c* k + 4-/- 

Eq. (40) gives 

. /f D 

^ - g 3 



(44) 



to ( C 5 hi + 3 hl+ sc,clh^c^)^A\ k, 



Thus 



Bo / ^ 1 j 2 - 3 s 

' = ^ & + 3C ' c > k +3C ' C ^ Ci 3 



(45) 



Inserting Eq. (45) into Eq . (44) 



4-4 4 + 6 * 4 + 6 / - ^ 



where 



SjL — b<> / c,/j 

3 C t ' c i ^ 3 bo 



a 



Q A 



-i>,C 



( 46 ) 
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■p — ;££/ £2 Q kjL s- A c i G 
^2 ' c e A 

8, = -SlJip . — + Ac s -b,C 
' c ( A 

In Eq. (46) all coefficients except b Q and can be 
obtained from hydrod-namic derivatives of a given ship. 
Calculator of b Q and b^ 

For 3 real negative roots of the characterisitic 
3 2 

equator, As + Bs + Cs + D = 0 , the point M on Mitrovics 
chart (see Ref. 7) should be in £ = 1.0 region 





For 3 repeated roots, M should be the cusp point. 
At the cusp point, 



J b 0 

d bj 

Tbi 



undefined 




Therefore 



dk°- = 0 W vT =0 

du)t 
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For the third order characteristic Eq . with 



1.0 



So, 



Note : 



j^o — U)t -2 UJ-fc 

b, - ZH>i -3 LOt 

- 2- bait = 0 

t 

<H U) f = -5T 

b, = 3 anc < k = If 

(see Ref. 7 for more details) 

It is difficult to find the exact cusp point by 
the digital computer due to the non-integer 
values of b_ and b-, . 
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SAMPLE 

MJM=NUM + 1 

CALL DRWG( 1 » 1 fNUMtYAXCtXAXCI 
CONTRL FINTIM= 29 .fDELT= 0 . 01 ?DELS =0 



PRINT . 50 ,XAXC, YAXC ,PSIC , VC »RC, TMCOEL 
TERMINAL 

CALL 6 NDRWCNPLOT) 

GO TO( 1 , 2 f 3 , 4 » 5 » 6 , 7 , 8 ) ,NPLGT 
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COMPUTER PROGRAM #3 
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